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We give a new separability criterion, a necessary condition for separability of N- 
partite quantum states. The criterion is based on the Bloch representation of a iV-partite 
<^ quantum state and makes use of multilinear algebra, in particular, the matrization of 

' tensors. Our criterion applies to arbitrary iV-partite quantum states in H = H dl <S> H d2 ® 

• • -<g)TC dN . The criterion can test whether a iV-partite state is entangled and can be applied 
to different partitions of the iV-partite system. We provide examples that show the ability 
of this criterion to detect entanglement. We show that this criterion can detect bound 
entangled states. We prove a sufficiency condition for separability of a 3-partite state, 
straightforwardly generalizable to the case N > 3, under certain condition. We also give 
a necessary and sufficient condition for separability of a class of iV-qubit states which 



includes iV-qubit PPT states. 



Keywords : Bloch representation of quantum state, separability criteria, matrization 
of tensors, PPT entangled states. 
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1 Introduction 



The question of quantifying entanglement of multipartite quantum states is fundamental 
to the whole field of quantum information and in general to the physics of multicomponent 
quantum systems. Whereas entanglement of pure bipartite states is well understood, the 
classification of mixed states according to the degree and character of their entanglement 
is still a matter of intense research [1,2]. A iV-partite state acting on TC = 7i dl <8> 'H d2 ® 
■ ■ • (g> 7i dN is separable [3] ( or fully separable) if it can be written as a convex sum of 
tensor products of subsystem states 

N 

P = ^^pL 1) ®pl 2) ---®pL iY) = ^^(g)p2' ) , P W >0; Y,P™ = 1 - (!) 

w w j=l w 
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A state is called k separable if we can write 



w 



,(Ofe) 



(2) 



w 



where a^; % = 1, 2, ... ,k are the disjoint subsets of {1,2,..., N} and acts on the 
tensor product space made up by the factors of H labeled by the members of a*. The 
understanding of multipartite entanglement has progressed by dealing with some special 
classes of states like the density operators supported on the symmetric subspace of 7i 
[4]. A lower bound on concurrence on the multipartite mixed states is obtained [5]. K. 
Chen and L. Wu have given a generalized partial transposition and realignment criterion 
to detect entanglement of a multipartite quantum state [6]. 

There are two definitions commonly used for the entanglement of multipartite quantum 
states, the one from Ref. [7] (ABLS) and the one introduced in [8] (DCT). In DCT, all 
possible partitions of N parties are considered and it is tested for each partition if the state 
is fully separable there or not. A state is called N partite entangled if it is not separable 
for any partition. If a state is separable for a bipartite partition, it is called biseparable. 
In ABLS, a state is called biseparable if it is a convex combination of biseparable states, 
possibly concerning different partitions. A V-partite entangled state is one which is not 
biseparable. 

In this paper we derive a necessary condition for the separability of multipartite quan- 
tum states for arbitrary finite dimensions of the subsystem Hilbert spaces and without 
any further restriction on them. The criterion is based on the Bloch representation of a 
multipartite quantum state, which has been used in previous works to characterize the 
separability of bipartite density matrix, in particular, our work is a generalization of de 
Vicente's work on bipartite systems [9]. We make use of the algebra of higher order 
tensors, in particular the matrization of a tensor [10,11,12,13,14,15,16,17]. 

The paper is organized as follows. In section II we present the Bloch representation of 
a iV-partite quantum state. In section III we obtain the main results on separability of a 
iV-partite quantum state. In section IV we give a sufficient condition for the separability 
of a 3-partite quantum state generalizable to the case N > 3. In section V we investigate 
our separability criterion for mixed states, in particular, bound entangled states. Finally 
we summarize in section VI. 



2 Bloch Representation of a TV-Partite Quantum 



Bloch representation [18,19,20,21,22] of a density operator acting on the Hilbert space of 
a d-level quantum system C d is given by [9] 



State 




(3) 
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Eq.(3) is the expansion of p in the Hilbert-Schmidt basis {Id,Xf,i = 1, 2, . . . , d 2 — 1} 
where Aj are the traceless hermitian generators of SU(d) satisfying Tr(Xi\j) = 25^- and 
are characterized by the structure constants of the corresponding Lie algebra, fijk,gijk 
which are, respectively, completely antisymmetric and completely symmetric. 



— ~^n^ d + ifijk^k + 9ijk^k 



s = (si, S2, ■ ■ • , Sd2-i) in Eq.(3) are the vectors in IR d2_1 , constrained by the positive 
semidefiniteness of p, called Bloch vectors [21]. The set of all Bloch vectors that consti- 
tute a density operator is known as the Bloch vector space ^(IR^ 1 ). The problem of 
determining BiR^- 1 ) where d > 3 is still open [19,20]. However, for pure states (p = p 2 ) 
the following relations hold. 



(4) 



S 2 



d(d-l) 



SiSjgijk = (d - 2)s k 



(5) 



D<2 2 -1 



where ||.||2 is the Euclidean norm in I 

It is known [23,24] that i^IR^ 2-1 ) is a subset of the ball D^R?' 1 ) of radius R 

yj ^ d( - d ~^ ; which is the minimum ball containing it, and that the ball D r (IR d2_1 ) of radius 
is included in ^(R^- 1 ). That is, 



r = 



2(d-l) 



D r 



l d2 - r ) C B(R d ' 2 ^ 1 ) C D R { 



»<i 2 -i 



(6) 



In order to give the Bloch representation of a density operator acting on the Hilbert 
space C dl (g>C d2 ®- • -®C dN of a iV-partite quantum system, we introduce following notation. 
We use k, ki (i — 1, 2, • • • ) to denote a subsystem chosen from N subsystems, so that k, 
ki (i — 1, 2, • • • ) take values in the set Af = {1, 2, • • • , N}. The variables a k or a kl for a 
given k or ki span the set of generators of SU(d k ) or SU(d ki ) group (Eqs.(3) and (4)) for 
the kth or kith subsystem, namely the set {Ai fc _, A2 fc ., • • • , A d 2 _ x } for the kith subsystem. 
For two subsystems k\ and ki we define 



= (J* ® J < 



= (J, 



<*2 



a£ 2) = ® J < 



<*2 



A a , <g> Jd t 



A afci <g>J, 



h N ) 



K k0 ® id fc , +1 ® • • • <8> h N ) 



K k2 ®h k2+1 ®h N ) (7) 



where A afci and A afc2 occur at the k\th and kith places (corresponding to fcith and A; 2 th 
subsystems respectively) in the tensor product and are the a^th and afc 2 th generators of 
SU (d kl ), SU (dk 2 ), ttfci = 1,2,..., — 1 and a k2 — 1, 2, . . . , d\ 2 — 1 respectively. Then 
we can write 
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p = wr^ 1 ^ + E E + E E v^t^ + • • • + 

fc fc feeA/" Ofe {fei,fc2} Q fei«/= 2 

V V t \(ki)\(k2) ...\(fcM), ... V f A (l) A (2)... A (JV)-j. 

{fci,fc2,-" ,feM> a fc 1 a fe 2 '" afe M aiora— ajv 

(8) 

where s( fc ) is a Bloch vector corresponding to fcth subsystem, s( fc ) = [s ak \a k =i which is 
a tensor of order one defined by 



S a k 



rr|>A«] = ^Tr[p fc A Qfe ], (9a) 



where p& is the reduced density matrix for the /cth subsystem. Here {ki, k%, ■ ■ ■ , Ivm}, 2 < 
M < N, is a subset of J\f and can be chosen in (^) ways, contributing terms in the 
sum J2{k! k 2 ■■■ fc M } * n E<4-(8), eacn containing a tensor of order M. The total number of 
terms in the Bloch representation of p is 2^. We denote the tensors occurring in the sum 

£ {fel , fe2 ,-,*M}> (2 < M < N) by r^-*».-. fc -> = [V^-^mI which are defined ^ 

_ 4l4 2 ■ ■ ■ d kM r Uk^U^) . . . A (fe M )l 

= dkld % M ' dkM Tr[ Pklk2 ...UK kl ® A afc2 ® • • • ® A Qfc J] (96) 

where p klk2 ... kM is the reduced density matrix for the subsystem {/ci/^ • • • k M }. We call 
The tensor in last term in Eq. (8) T (7V) . 



3 Separability Conditions 

Before we obtain the main results we need following definition. Throughout the paper, 
we use the bold letter for vector and normal letter for components of a vector, matrix and 
tensor elements. 

A rank-1 tensor is a tensor that consists of the outer product of a number of vectors. 
For Mth order tensor and M vectors this means that = 

uf^uf 2 . . . for all values of the indices. This is concisely written as = o 

u ( 2 )o-.. ou W [17,11]. 

Also, given two tensors and of order M and A" respectively, with dimensions 
7x x I 2 x • • • x I M and J\ x J 2 x • • • x J N respectively, their outer product is defined as 
[16,10] 

(T 1 - ' O 5' ^)iii2—iMjih—3N ~ thi2-iM S ji32-jN (^) 
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Proposition 1 : A pure iV-partite quantum state with Bloch representation (8) is 
fully separable (product state) if and only if 

<J-{kl,k2,- ,k M } = s (fci) s (k2) o . . . o s (feM) ^11^ 

for 2 < M < N. In particular T (7V) = s« o S ( 2 ) o • • • o s^ holds. Here {h, fc 2 , . . . , C 
{1,2,..., iV}, and s( fc ) is the Bloch vector of kih subsystem reduced density matrix. 

Proof : Notice that Eq.(8) can be rewritten as 



{ki,k 2 } a k 1 a k 2 

+ •••+ [^aicw-ajv ~~ S Qi S Q2 • • • s Qjv]-^Li ' ' ' ^aj}- (^) 

For full separability, the sum of all the terms apart from the first term must vanish. 
Note that for every subsystem k — 1, 2, . . . , TV the set {Id, Xf, % — 1, 2, . . . , d\ — 1} forms 
an orthonormal Hilbert-Schmidt basis for the fcth subsystem. Hence A^; A^A^ . . . ; 

A^A^ • • • Xa/^; • • • ; Aa^A^ • • • A^ are the vectors belonging to the orthonormal prod- 
uct basis of the Hilbert-Schmidt space of the whole iV-partite system. By orthonormality 
of the tensor product of A's occurring in different terms, the required sum will vanish if 
and only if coefficients of each term vanish separately, that is if and only if 

toc kl a k2 -a kM = S aki S ak2 . . . Sa kM ', 2 < M < N, 

that is, 

j-{ki,k2,- ,k M } — s (fci) o s (fo) o • • • o s( fciV/ ) ; 2 < M < N. □ 

In fact, the condition (11) for all N parts is enough to decide the separability of pure 
iV-partite quantum states, as the following proposition shows. 

Proposition la : A pure TV-partite quantum state with Bloch representation (8) is 
fully separable (product state) if and only if 

where s( fc ) is the Bloch vector of kth subsystem reduced density matrix. 
Proof : Suppose p is a product state p = p 1 ® p 2 ® • • • <8> Pn- Then, 

G?lG?2 • • • diN rn r/ \/\ \ M 

t ai a 2 ...a N = Tr i\P^ ® P2 ® • • • ® pjv)(A Ql <g> A„ 2 <g> • • • <g> A Qjv )J 



= dld2 2N ' dN TrKpiKj ® (p 2 A Q2 ) <g> • • • <g> (pjvA ajv )] 
= dld2 2N - dN [Tr( Pl X ai )Tr(p 2 \ a2 ) • • • Tr{p N K N )} 

Suppose the condition holds, that is, [s^ o o • • • o s^] aia2 ... aN = t ai0l2 ... aN . Then, 

[s {1) o s (2) o • • • o s (Ar) ] QlQ2 ... ajv = dld2 2N ' dN [Tr(p 1 X ai )Tr(p 2 X a2 ) ■ ■ -Tr(p N X aN )] 

d\d 2 ■ ■ ■ d]\f .. . . . \ \i 

— Tr [(piA Ql ) <g> (p 2 A Q2 ) ® • • • ® [pNX aN )\ 



2 N 

d\d 2 ■ ■ ■ dj\f 



2 N 



Tr[(pi ® p 2 <g) • • • <8> Pat) (A Ql <g> A Q2 <g> • • • <g> A 



d\d 2 . . . d]\f r . \ \i 

= t ai0l2 ... aN = ^ Tr[p(X ai <g> A Q2 <g> • • • <g> A ajv )J. 

The equality 

7V[(pi <g> p 2 ® • • • ® pjv)(A ai <g> A a2 <g> • • • <g> A ajv )] = Tr[p(A ai <g> A a2 ® • • • <g> A Q J] 

is satisfied for all elements in the ortho normal basis {^>k=i^a k } , < < g? 2 — 1, (a k = 
for Jd fc ) where {X ak } are the d\ — \ generators of SU(dk). This means that the joint 
probabilities obtained from the ensemble of measurements of (A ai •••A QJV ) for states p 
and p = pi ® p 2 ® • • • ® Pn are equal. This implies 

p = pi <8> P2 <E> • • • <E> Pn- 

□ 

Note that this criterion is easily amenable with experiments. In order to check it for 
an element of we have to measure the corresponding generators on each subsystem 
and then check whether the product of the averages equals the average of the products. 

Thus in order to check whether a given pure state is a product state we have to check 
whether = os' 2 ' o • • •os^, where the Bloch vectors s^, s^ 2 \ . . . , can be con- 
structed from the reduced density matrices p±, p 2 , ■ ■ ■ , pn for subsystems 1, 2, • • • ,N (s ak = 
f Tr(p k X ak ), k e {1, 2, • • • , N}, see Eq.(9a)). 

In the case of mixed states we can characterize separability from the Bloch represen- 
tation point of view as follows. 
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A N -partite quantum state with Block representation (8) is fully separable if and only 
if there exist vectors uffl G IR d fe -1 satisfying Eq.(5), and weights p w satisfying < p w < 1 
and YliwPw — 1 su °h that 

T^ = Y. Pw OL^\ s<*> = X> w t*S> (13a) 

w w 

and 

R 

T {kite,-,kM} = J2p w Ofii (136) 

w 

for 2 < M < N ; for all subsets {k 1: fc 2 , . . . , k M } C {1,2,..., N}, 

where s( fe ) is the Bloch vector of the mixed state density matrix for kth subsystem and 
uffi represent the Bloch vector of the pure state of the fcth subsystem contributing to the 
wth term in Eq. (1). 

This follows from proposition 1 and Eq. (1). However, in view of proposition (la), the 
necessary and sufficient condition is given by Eq.(13a), so that Eq.(13b) can be dropped. 
The above result can not be used directly, as it amounts to rewriting Werner's definition 
of separability in a different way. However, it allows us to derive a necessary condition 
for separability for iV-partite quantum states. 

We need some concepts in multilinear algebra. Consider a tensor G ^hxi 2 x-xi N 
, where I k — d\ — 1. The nth matrix unfolding of T^ N "> (n = 1, 2, • • • , N) [10] is a matrix 
G ^nx(i n+1 i n+2 ...i N hh-in-i)^ T W con t a i ns the element t hi2 ... iN at the position with 
row index i n (i n — 1, 2, • • • , I n ) and column index 



(in+l ~ l)^n+2-^n+3 • • • -TtV-Tl^ • • • In-1 + (^n+2 — l)^n+3-^n+4 • • • ■ ■ ■ In-1 

H \- (in - l)hh ■ ■ ■ I n -i + (h - ±)hh ■ ■ ■ h-i + {k - ^-)hh ■ ■ ■ h-i H h i n -i- 

For n — 1, we take the last term i n _i = i = i^. This ordering is called backward cyclic 
[16]. To facilitate understanding, put N points on a circle and label them successively 
by zi , z 2 , • • • ,In- The consecutive terms in the expression for the column index in T^y* 
corresponding to t^,^,... ,i N become quite apparent using this circle. 

For G R /lX/2X/ 3 the matrix unfolding T^) contains the elements t iliai3 (i^ = 

1, 2, • • • , Ik] k — 1, 2, 3) at the position with row number %\ and column number equal to 
(i2 — 1)13 + 23, T}J contains t^i^ at the position with row number 12 and column number 
equal to (i 3 — + i\ and T^) contains t ili2i3 at the position with row number i 3 and 
column number equal to {i\ — l)I 2 + ii- 
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^123 = ^322 = 0. The matrix unfolding T^) is given by 



T (3) 




As an example [25], define a tensor T (3) G M 3x2x3 , by t ul = t 112 = t 2 u = -t 2 i2 = 1, 

^213 — ^311 — ^313 — ^121 — t\22 = ^221 = — ^222 = 2, £223 = ^321 = ^323 = 4, tn3 = £312 = 

,(3) 

(1) 

/ 1 10 
1 -1 2 
\ 2 2 

Note that there are N possible matrix unfoldings of T^ N \ The matrix unfolding is 
called the matrization of the tensor [10,17]. We can now define the Ky Fan norm of the 
tensor 7~( N ) (of order N) over N matrix unfoldings of a tensor, as 

\\T {N) Wkf = max{ \ \T$ \ \ KF }, n = 1, . . . , N; (14) 

where | \T^V \\ KF is the Ky Fan norm of matrix T^V defined as the sum of singular values 

of [26]. It is straightforward to check that \\T (m \\ kf defined in (14) satisfies all the 
conditions of a norm and is also unitarily invariant [9,26]. 

The tensors in Eq.(13a) are called Kruskal tensors with the restriction < p w < 
1, ^2 w p w = 1 [14,16]. We are interested in finding the matrix unfoldings and Ky Fan 
norms of occurring in Eq.(13a). The kth matrix unfolding for Kruskal tensor is [17] 

T (A0 = [/ W E ( [/ (AT) jjiN-l) . . . ^(fc+1) ^(fc-1) . . . j^T (15) 

Here = [u[ k) u { 2 k) . . . uf] G M /fcXR ; k = 1, 2, ... N and R is the rank of Kruskal 
tensor [14,12,17], i.e. the number of terms in Eq.(13a). uf^ is a vector in M /fe and is the 
ith column vector in the matrix U^ k \ X is the Rx R diagonal matrix, X =diag[pi . . -Pr\. 
The symbol denotes the Khatri-Rao product of matrices [17] U G R /xi? and V G M JxR 
defined as[/"07= [u^V! u 2 ®v 2 ... Ur® v R ] G R IJxR where U; and Vj, % — 1, 2, . . . R 
are column vectors of matrices U and V respectively. Eq.(15) can be rewritten as 

T$ = U^[v?v? ■ ■ • = U^YV^ (16) 

where v-^; % = 1, 2, . . . , R are the column vectors of the matrix 

V {W e m I N I N -ilN- 2 ...I k+1 I k - l ...hxR and v f ) = u f) U (N-1) u (N-2) . . . u (fe+l) 

(k) 

^=Ep»uff)vJ) T ; A; = 1,2,..., AT. (17) 

UI = 1 



u[ k ^ ® • • • <g> uf \ Using Eq.(16) we can write Tffi as 



Theorem 1 : If a iV-partite quantum state of d\d 2 ■ ■ ■ dm dimension with Bloch 
representation (8) is fully separable, then 



\\T (N) \\KF<J^ =1 d k (dk-l). (18) 
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Proof : If the state p is separable then T( N ) has to admit a decomposition of the 
form Eqs.(13) with ||u^' ) ||2= \J dk ^~^ -, k = 1,2, ...,N. From definition of KF norm of 
tensors, Eq.(14), 

\\T (N) \\kf = max{ \ \T^ \ \ KF } ; k = 1, . . . , N. 

From Eq.(17), 

\\T {N) \\kf = max{ \ | £ Pw u ( w k ^ T \ \ KF } ; k= 1, . . . , N 



< max{^ Pw \\ii { ^ )T \\ KF } = max{J2p w y^kdk(dk - l)\\u^v^ T \\ KF } 

w w 

where ul fe) , vi fe) are unit vectors in M^ _1 and M^" 1 <g> <g> • • • <g) ]R d fc+i _1 ® 

M^Li^ 1 ® . . . <g) R^?- 1 respectively, so t hat Wu^^ \\ KF = 1 for all k = 1, 2, ... , TV. 
Using = 1 we get | |TW 1 1 ^ < ^^4(4-1). □ 

For a subsystem we get, 

Corollary 1 : If the reduced density matrix of a subsystem consisting of M out of 
N parts is separable then \\T^ k ^- ^\\ KF < \J ^Uf =l d k {d k - 1). 

The negation of the above condition, that is, > yj ^wll k=1 d k (d k — 1), is a 

sufficient condition of entanglement of iV-partite quantum state. This leads to a hierarchy 
of inseparability conditions which test entanglement in all the subsystems. For N = 2 

the condition ||T^||xf < \J ^4(4 — 1)4(4 _ 1) has been shown in Ref. [9], to be a 
sufficient condition for entanglement associated with any bipartite density matrix. Note 
that for iV-qubits, 4 = 2, % — 1, 2, . . . , TV, the above criterion becomes, for a separable 
state, IIT^IIa-f < 1. 

Consider a TV qudit system 7i s = ®^ =1 7i k in a state p, supported in the symmetric sub- 
space of 7i s . It is straightforward to see that all the tensors in the Bloch representation of p 
are supersymmetric, that is (see Eqs.(8) and (9)), t akiCXk2 ... akM = fp( afcl )p(a fca )-p(a fcM ), 2 < 
M < N, where P is any permutation over indices {a kl , a k2 , • • • , ot kM }. We have, neglect- 
ing the constant multipliers, 



t akl a k2 -a kM ~ Tr[p klk2 ... kM \ aki <g) <g) • • • <g) A afcj J 

= Tr[p klk2 ... kM PP T X aki <g) A afc2 ® • • • <g> A QfcM PP T ] 

= Tr[(P T p fclfc2 ... fcM P)(P r A afci ® A Qfc2 ® • • • ® A QfcM P)] 

= Tr[p klk2 ... kM (X P{aki) <g> Ap Ka ) (8) ••• (8) A PKm) )] 

= ^P(a fcl )PK 2 )-P(a fcM ) 
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where P is the appropriate permutation matrix permuting the A matrices within the 
tensor product [26], P T being the transpose of P satisfying P T = P^ 1 . In particular 
is supersymmetric. All matrix unfoldings of a supersymmetric tensor have the same set 
of singular values [10] and hence the same KF norm. Thus, for a iV-qudit system in a 
state supported in the symmetric subspace, it is enough to calculate the KF norm for any 
one of the N matrix unfoldings to get max{||T^||x.F}- 

4 A Sufficient Condition for Separability of a 3- 
Partite Quantum State 

Consider the Bloch representation of a tripartite state p acting on H dl <8> 7~i d2 <8> H ds , 
di < d 2 < d 3 . 

Oil Ct2 &3 aiOi2 

+ V« A (1) A (3) + Vt A (2) A (3) + A (1) A (2) A {3) (19a) 

010:3 0203 010203 

where r, s and q are the Bloch vectors of three subsystems respectively , T^' u ^ = 
[t a av ] the correlation matrix between the subsystems /x, v\ v\ C {1,2,3} and = 
[fo lCt2Q3 ] the correlation tensor among three subsystems. Before stating proposition 2, we 
need the following definition and result. 

Kruskal decomposition of a tensor 

j=i 

is called completely orthogonal if (v%\ uf^) = 5m, i — 1, 2, • • • , N; k, I — 1, 2, • • • , R 
[13], where (, ) denotes the scalar product of two vectors. If has completely orthog- 
onal Kruskal decomposition, then it is straightforward to show that 

11^11^ = X>, (20) 
j'=i 

where R is the rank of and £j, j = 1, 2, • • • , R are the coefficients occurring in the 
completely orthogonal Kruskal decomposition of T^ N \ In the proof of proposition 2, we 
assume that completely orthogonal Kruskal decomposition of T^ k \ k > 2 is available. A 
completely orthogonal Kruskal decomposition may not be available for an arbitrary tensor 
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[13]. The general conditions under which the completely orthogonal Kruskal decomposi- 
tion is possible is an open problem. We conjecture that completely orthogonal kruskal 
decomposition is available for all tensors in the Bloch representation of a quantum state, 
but we do not have a proof. As it stands, this issue has to be settled case by case. 

Proposition 2 : If a tripartite state p acting on H dl ® H d2 <8> H d3 , di < d 2 < 
d 3 , with Bloch representation (19a), where has the completely orthogonal Kruskal 
decomposition, satisfies 



\KF 



8 ( dl _ i)( d2 _ i)( d3 _ i) 
+ ^ l|T|1 ^ " h (21) 



then p is separable. 



Proof : The idea of the proof is as follows. 



(i) We first decompose all the tensors in the Bloch representation of p as the completely 
orthogonal Kruskal decomposition in terms of the outer products of the vectors in the 
Bloch spaces of the subsystems (coherence vectors). 

(ii) We prove that we can decompose p using the Kruskal decompositions described 
in (i) above, as the linear combination of separable density matrices, which is a convex 
combination if the coefficient of identity is positive. This condition is the same as the 
condition stated in the proposition. 

Let T^'"}; {//, u} C {1,2,3} in Eq.(19a) have singular value decomposition = 
E^a^a^f ; with ||a?°|| 2 = ||af } || 2 = 1 , for {/i, z/} C {1,2,3} and let T in 
Eq. (19a) have the completely orthogonal Kruskal decomposition T = ^Uj o Vj o wj 
[17,14,27] with ||uj|| 2 = ||Vj|| 2 = ||Wj|| 2 = 1. We define 

^ = 4^^ , {1,2,3} 
so that we can rewrite 

= ^Z^zA ^a^f (22a) 

Similarly, we define 
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G i = \lwh) U J ! = ^/ippi) v i ' = ^/a^TJ w i> so that we can write 



y ^ ^ Vj Wj (226) 

v j 

If we substitute Eqs.(22a) and (22b) in p Eq.(19a), we get 



+ 



p = d^d 3 ^ 3k=iIdk + E r ^ A S } + E s °* A S + E SosAg 

ai Q!2 <*3 

E E J 4(d " "ffi " 1} £^i, M w(i, M UAgAg 
J a(dl = ^j^^ E E««<).ft)-(*i)-W»S d») 

V "' " ai (was 7 



+ 



The coherence vectors a- occur in occur m Uj occur m 

D r (R^ _1 ), Vj occur in /^(R^ 1 ) and w,- occur in D r (R d ^ r ) (see Eq.(6)), so that they 
correspond to Bloch vectors. 

We can decompose p Eq.(19b) as the following convex combination of the density 
matrices 

Pj > P'i > Pj > Pj" i ft > ^ > T i > T i > ^ > ^ i Pr > Ps > Pg aIld d^k Id l<kdz\ 



n V" - 1)(^2 - l)(rf 3 -1) ^^ ,„/ , // | ///N | \- / 4(^l-l)(4-"Ty ^^ , ,x 

p = E y ^ j(pj+Pj+Pj+Pj y — ^ — ¥ (ft+^) 

+ Ey — ^ — y^+^+Ey — ^ — y^+tj 



'8(dl-l)(d 2 -l)(d 3 -l) ||/T . 1| N / dl d 2 d3 



did 2 d 3 did 2 dz 
where pj in Bloch representation is 



r\\K F )-P^-. (23) 
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rf, 



jj-r{ ®Li ^ + £(fl,)«A<'> + £(v,)„Ag,' + £(w,)«,A£ 

CK1 CK2 Q3 

0102 0103 020:3 

^ (u J 0« 1 (v i ) oa (w i )«,A£)AgAg) 



+ 

010203 

= jhk {Idl + ^ (G ^ ) - A - ) ^ + E^Ag) ® (J* + £(*;)«,Ag). (24) 

Ql O2 O3 

Note that ||T||xf in Eq.(23) is denned via Eq.(20), which is based on completely 
orthogonal Kruskal decomposition of T . 

The Bloch vectors, correlation matrices and correlation tensors of the density matrices 
Pj , Pj , Pj , ft" ; Qi, Qi , n , t- , TTi , 7T- ; pr , ps , p q are 

For Pj , 

~ rp{l,2} ~ ~t t-.{1,3} - ~T t^{2,3} - ~t 

For pj, 

r', = Uj , s', = -v, , q', = -w, , T'^ = -ujvj , if ' 3} = -u,wj 
lf 2 ' 3} = Vjwj , T'j — Uj o vj o Wj. 

For p" 



r 3 — U 3 ' J J ) Q. j W J ) J j - U ;? V j > J j - U J W j 

rp»{2,3} _ _y^J , T"j — Uj O Vj O Wj. 

For pj, 

T-,///{2,3} ~ ~ r n-ill 

For Qi, 

r e _ 3 (i) _<? _ 3 (2) _e _ n T c{i,2} _ 5 (i) 3 (2)T T <?{1,3} _ n 
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T 8{2,3} = Q r e = 
For g'i, 

r e' _ =(1) Jj _ =(2) o' _ n T Q'{1,2) _ ~(1)~(2)T 
^'{1,3} = q ) ^'{2,3} = Q %e ' = Q 

For rw, 



r[ = af } , s[ = , q[ = af , if 1 ' 2 * = , 7? 

^{2,3} = q ? ^ T = q 

For r/ 



{1,3} _ S (1) 5 (3)T 



r r' _ =(1) „r' _ n _r' _ s(3) ^'{1,2} _ Q ^'{1,3} _ 3 (1) 5 (3)T 

r i — _ d i 5 S i — U ) — _ d i > J j — U J J i — d j d i 

T r'{2,3} 

For 7r, 



T[' {2 ' 3} = =0. 



if = , = af } , q? = af } , T^l, 2} t = , 7f {1 > 3 > = 

^{2,3} = ~(2)~(3)T iVt = Qm 

For 7r', 

if = , sf = -af ) , qf = -af > , 7f ^ = , if™ = 

^'{2,3} = 5 (2)-(3)T ^ v > = q 

For p r , 



r r = ^2(fc)Mh , s r = , q r = , T r { ^ } = ; V{p, z/} C {1,2,3} , % = 0. 
For p s , 



r * = > s * = a/27^™ > q S = , T s { ^ } = ; V{/i,z/} C {1,2,3} , % = 0. 



2(da-l) ||s|| 2 

For p ? , 

r g = 0, s, = 0, q, = V / 5 Z TJ^, = ; V{p, z/} C {1,2,3}, T s = 0. 

If we write all matrices p'j , p" , p"' ; Qi , ^ , Tj , , 7Tj , 7r^ ; p r , p s , p g (as we have 
done for pj in Eq.(24)) in the Bloch representation and substitute them in Eq.(23) we get 
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p as in Eq.(19b). 

To understand this let us see how the first term in Eq.(23) adds up to give the 
last term in Eq.(19b). The definition of pj, p'p p", p'" (denoting the Bloch vectors 
by si, S2, S3, S4, ....) can be summarized in the tabular form 



Table 1 

Correspondence between the first term in Eq.(23) and the last term in Eq. (19b). 





Si 


S 2 


S3 


S1S2 


S1S3 


S2S3 


S1S2S3 


Pi 


Uj 


V 3 


Wj 


UjVj 


UjWj 


VjWj 


UjVjWj 


p'l 


Uj 




-Wj 


-UjVj 


-UjWj 


VjWj 


UjVjWj 


p" 


-Uj 




-Wj 


-UjVj 


UjWj 


-VjWj 


UjVjWj 


p'i 


-Uj 




Wj 


UjVj 


-UjWj 


-VjWj 


UjVjWj 



The contribution of each column to pj + p 3 + p" + p"' is zero except the last column 
which reproduces the last term in Eq.(19b). We can get the contributions of each term 
in pj, p'j, p", p"' to their sum by just keeping track of their signs. Thus we only need the 
following table (dropping j) 



Table 2 

Contributions of various terms in p,p',p",p"' to their sum. 





Si 


S2 


S3 


SlS 2 


Si S3 


S2S3 


S1S2S3 


p 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


p' 


+ 










+ 


+ 


p" 




+ 






+ 




+ 


p'" 






+ 


+ 






+ 



In the same way, the contributions of the terms involving g, r, ir are obtained by using 
the table corresponding to table 2 for the bipartite case [9]. g, r, it which contain tensors 
of order two correspond to three 2-partite subsystems 12,13 and 23 . The corresponding 
tables are 



Table 3 

Contributions to g + g' 
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S 2 


S3 


S1S2 


Si S3 


s 2 s 3 


S1S2S3 


g 


+ 


+ 





+ 











d 









+ 












Table 4 

Contributions to r + t' 





Si 


s 2 


s 3 


sis 2 


S1S3 


S2S3 


S1S2S3 


T 


+ 





+ 





+ 








t' 












+ 









Table 5 

Contributions to n + it' 





Si 


s 2 


s 3 


sis 2 


S1S3 


S2S3 


S1S2S3 


7T 





+ 


+ 








+ 





7T' 















+ 






Tables 2, 3, 4, 5 encode the procedure to construct the possible separable state given 
in Eq.(23). 

We now note the following points 

(i) If the condition (21) holds, then the coefficient of the matrix Id^dz in 
Eq.(23) is positive which ensures that the decomposition (23) of p is posi- 
tive semidefinite. 

(ii) By virtue of Eq.(6), all the coherence vectors occurring in p'j , p" , p"- ; 
Qi , Qi , Ti , t[ , 7Tj , 7T- ; p r , p s , p q belong to the corresponding Bloch 
spaces. 

By (i) and (ii) we conclude that p'j , p'- , p"- ; q { , q\ , n , t[ , 7Tj , 7r- ; p r , p s , p q 
constitute density matrices. Further, all these matrices satisfy condition (11) so that, via 
proposition 1, all these matrices correspond to pure separable states, equal to the tensor 
products of their reductions. Therefore, they constitute density matrices and they are 
separable and so must be p. □ 

We can generalize proposition 2 to the iV-partite case by constructing the tables 
successively for iV = 4, 5, 6, • • • . First note that the number of p s in the first term of 
Eq.(23) lifted to the iV-partite case is 2 N ~ 1 . For iV = 4 we have eight. The corresponding 
table is 
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Table 6 

Generalization of Table 1 to N = 4. 





Si 


S2 


S3 


S 4 


SlS 2 


S1S3 


S1S4 


S2S3 


S2S4 


S3S4 


S1S2S3 


S1S2S4 


p(l) 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


p (2) 


+ 


+ 






+ 










+ 






p (3) 


+ 




+ 






+ 






+ 






+ 


p (4) 


+ 






+ 






+ 


+ 






+ 




p (5) 




+ 


+ 








+ 


+ 








+ 


p (6) 




+ 




+ 




+ 






+ 




+ 




p (7) 






+ 


+ 


+ 










+ 


+ 


+ 


p (8) 










+ 


+ 


+ 


+ 


+ 


+ 







(Table 6. Continued) 



S1S3S4 


S2S3S4 


S1S2S3S4 


+ 


+ 


+ 


+ 


+ 


+ 




+ 


+ 




+ 


+ 


+ 




+ 


+ 




+ 






+ 






+ 



We see that the contribution of each column to the sum is zero except the last 

one corresponding to the Kruskal decomposition of occurring in the Bloch represen- 
tation of the given state p. For general case of iV-partite state we construct the table 
for pW, % — 1,2, • • • ,2 N ~ 1 as follows. First column consists of 2 N ~ 2 plus signs followed 
by 2 N ~ 2 minus signs. Second column comprises alternating 2 JV ~ 3 plus and minus signs. 
Continuing in this way upto 2 N ~ N = 1 we get alternating plus and minus signs in the 
(N — l)th column. We set the iVth column to ensure that there are zero or even number 
of minus signs in each row. Rest of the columns can be constructed by appropriate mul- 
tiplications. This procedure can be checked on table 6. We denote the sequence of such 
tables for TV = 2, 3, 4, • • • as T,, i = 2, 3, 4, • • • . 

The tables corresponding to (N — 1), (N — 2), ...,2 partite subsystems giving rise to 
the remaining terms in the equation (23), lifted to iV-partite obtained from 

T/v-i,Tjv_ 2 , ...,T 3 ,T 2 , exactly as described in the proof of proposition 2. In this way we 
can lift eq.(23) to the iV-partite case, with the total numbers of terms J^u^q 1 (^)2 N ~ 1 ~ l + l. 
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Once this is done, the rest of the proof for iV-partite case follows as in proposition 2. Thus 
we have 



Proposition (2a): If a iV-partite state p acting on H = H dl ®7i d2 <g> • • -®H dN , d x < 
d 2 < ■ • ■ < (In with Bloch representation (8), where all T^ k \ k > 2 have the completely 
orthogonal Kruskal decomposition, satisfy 



k V ak {fl , u} V a » av 

+ £ J 8(d "" 1) f- 1)(d ""^l|r'-'|U f + ... + 

(21) 

then p is separable. □ 
For a iV-qubit system Theorem 1 and proposition 2a together imply 
Corollary 2 : Let a iV-qubit state have a Bloch representation 



P 2? 



and let the tensor in the second term have the completely orthogonal Kruskal decompo- 
sition. Then p is separable if and only if ||T^||^ < 1. □ 



5 Examples 

We now investigate our separability criterion (18) for mixed states. We consider iV-qubit 
state 

p ( noly = l^rl+PMW, 0< P <1 (25) 

where \ip) is a iV-qubit W state or GHZ state. We test for N = 3,4,5 and 6 qubits. 
We get, 
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Table 7 

The values of p above which the states are entangled. 



\GHZ) 


\W) 


N 


p> 


p> 




0.35355 


0.3068 


3 


0.2 


0.3018 


4 


0.17675 


0.30225 


5 


0.1112 


0.3045 


6 



Entanglement in various partitions of W noisy state Eq.(25) is obtained by using 
(N — n) qubit reduced W noisy state 

P%£\w) = ^In-u + ^p\0 N - n )(0 N - n \ + ^^p\W N - n ){W N . n \ (26) 

For N = 6 and n = 2 we found that the state is entangled for 0.491 < p < 1. 
For N qutrits (d = 3) we test for 

P ( noly =1 -^r J + PMM (27) 

where = ^= Ylt=i • • •) ^ s ^ ne maximally entangled state for N qutrits. 
For N = 3 and N = 4 (qutrits) the state p£2 sy m -^1- (27) is entangled for 

0.2285 <p<l, N = 3 

0.2162 <p<l, N = A (28) 

The state 

Pno lsy = + p\i>){i>\ (29) 

where = |(|H2) + |123) + \2U) + |234>) in the space C 2 ® C 3 <g> C 4 is found to be 
entangled for 0.24152 < p < 1. 

All of the above examples involve NPT states. Now we apply our criterion to PPT 
entangled states for which PPT criterion is not available. 

We apply our criterion to the three qutrit bound entangled state considered by L. 
Clarisse and P. Wocjan [27], given by p c <g> \ip){ip\ where p c is the chess-board state given 
in [27] and \ip) is an uncorrelated ancilla. Our criterion detects the entanglement of this 
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state as ||T( 12 )|| = 3.75 > 3. Further, the four qutrit state p = (1 - 0)p c <g> p c + /3J/81 
considered by the same authors yields entanglement for < (3 < 0.2, after tracing out 
either subsystems 1 and 2 or subsystems 3 and 4. 

Now we consider the important example of the Smolin state [28,29], which is a four 
qubit bound entangled state given by 

i 4 

PTbcI = 4 E I^ASX^I ® IfeXfel (30) 
i=l 

where I^ab) anci IV'cd) are the Bell states. Pabcd ^as the Bloch representation 
ptmtoc* = ^(J® 4 + J2 3 i=1 af A ) so that Corollary 2 applies (note that the requirement of 
completely orthogonal Kruskal decomposition is trivially satisfied). We find for this state 
||^"^||a"f = 3 > 1 confirming its entanglement. 

Our last example is the four qubit bound entangled state due to W. Diir [30,31] 

tf E = \(\m\ + \ 1 E(Pi+Pi)) 
i=i 

where is a 4-party (GHZ) state , Pi is the projector onto the state which is a 
product state equal to |1) for party % and |0) for the rest , and Pi is obtained from Pi by 
replacing all zeros by ones and vice versa. We get HT^H^ = 1.4 > 1 confirming the 
entanglement of this state. 



6 Summary 

In conclusion we have presented a new criterion for separability of N partite quantum 
states based on the Bloch representation of states. This criterion is quite general, as it 
applies to all iV-partite quantum states living in 7i = 7i dl <8> 7i d2 <8> • • • <8> 7i dN , where, in 
general, the Hilbert space dimensions of various parts are not equal. Most of the previous 
such criteria had restricted domain of applicability like the states supported on symmet- 
ric subspace [4] or, are, in general, restricted to bipartite case. In proposition 2 we have 
given a sufficient condition for the separability of a tripartite state under the condition 
that the tensors occurring in the Bloch representation of the state have completely or- 
thogonal Kruskal decomposition. This result can be generalized to the iV-partite case. 
Via corollary 2 we give a necessary an sufficient condition to test the separability of a class 
of iV-qubit states which includes iV-qubit PPT states. Smolin state (30) is an important 
example in this class. The key idea in our work is the matrization of multidimensional 
tensors, in particular, Kruskal matrization. We have defined a new tensor norm as the 
maximum of the KF norms of all the matrix unfoldings of a tensor, which is easily com- 
puted. We have also shown that this norm can be calculated even more efficiently for 
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a iV-qudit state supported in the symmetric subspace. It will be interesting to seek a 
relation of this tensor norm with other entanglement measures. Again, the entanglement 
measures like concurrence known so far are successfully applied to pure states, bipartite 
or multipartite, while our tensor norm can be easily computed for arbitrary iV-partite 
quantum state. Finally, our result on full separability (proposition 1) of iV-partite pure 
states can be easily moulded for the fc-separability of an iV-partite pure state. In fact it 
is straightforward to construct an algorithm giving the complete factorization of the N- 
partite pure state (see the paragraph following the proof of proposition 1). It is also easy 
to see that theorem 1 can be applied to any partition of a iV-partite system via the Bloch 
representation in terms of the generators of the appropriate SU groups. Most important 
is the observation that all the tensors in the Bloch representation can be computed using 
the measured values of the basis operators {\ ak } so that our detectiblity criterion is ex- 
perimentally implementable. 
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